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ABSTRACT 


We present extensive calculations of radiative transition rates and electron 
impact collision strengths for Fe II. The data sets involve 52 levels from the 3d 
3d^4s, and conhgurations. Computations of A-values are carried out with 

a combination of state-of-the-art multiconhguration approaches, namely the rela¬ 
tivistic Hartree-Fock, Thomas-Fermi-Dirac potential, and Dirac-Fock methods; 
while the i?-matrix plus intermediate coupling frame transformation, Breit-Pauli 
R-matrix and Dirac i?-matrix packages are used to obtain collision strengths. We 
examine the advantages and shortcomings of each of these methods, and estimate 
rate uncertainties from the resulting data dispersion. We proceed to construct 
excitation balance spectral models, and compare the predictions from each data 
set with observed spectra from various astronomical objects. We are thus able to 
establish benchmarks in the spectral modeling of [Fe II] emission in the IR and 
optical regions as well as in the UV Fe II absorption spectra. Finally, we provide 
diagnostic line ratios and line emissivities for emission spectroscopy as well as 
column densities for absorption spectroscopy. All atomic data and models are 
available online and through the AtomPy atomic data curation environment. 

Subject headings: atomic data - quasars: absorption lines - quasars: individual 
(QSO 2359-1241) - ISM: individual objects (HH 202) - ISM: individual objects 
(Orion) - ISM: individual objects (ESO-Ha 574, Par-Lup 3-4) - ISM: jets and outflows 
- ISM: lines and bands - stars: pre-main sequence 
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1. Introduction 


Reliable quantitative spectral modeling of singly ionized iron (Fe II) is of paramount 
astrophysical importance since various fundamental research lines—e.g. active galactic 
nuclei (AGN), cosmological supernova light curves, solar and late-type-star atmospheres, 
and gamma-ray-burst (GRB) afterglows—depend on such models. This ion has gained even 
more attention in recent years with the advent of large-scale observational surveys as well 
as deeper and high-resolution spectroscopy. 


Fe II spectral modeling hrstly requires the detailed treatment of electron impact 
excitation of metastable levels followed by spontaneous decay through dipole forbidden 
transitions. The computation of accurate electron impact collision strengths and A-values 
has proven to be cumbersome despite many efforts over several decades. The difficulty 
in describing the Fe II system arises from the complexity of the effective potential acting 
on the 4s, 3d, and 4p electrons. In practice, the wave function of the atomic system is 
approximated by an anti-symmetrized product of one-electron radial functions determined 
from the effective potential 


= V{nlr) -|- 


/(/ + !) 
2r2 ^ 


( 1 ) 


where V{nlf) is the electrostatic potential arising from the nucleus and the {N — 1) electrons 
of the ion, the second term being the centrifugal energy for an electron with orbital angular 
momentum quantum number /. For an effective potential with asymptotic form —2/r the 
centrifugal term takes the for m of a positive barrier for I > 2, and the effective potential 


becomes a two-well potential (iKaraziva 


19811). The two potential wells are very different 


from each other, the inner well is determined by many-electron effects while the outer is 


^At present at Astrophysique et Spectrocopie, Universite de Mons - UMONS B-7000 


Mons, Belgium 
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mostly hydrogenic. Therefore, slight variations in the potential morphology can lead to 
large changes in electron localization. For this reason, the atomic structure is very sensitive 
to orbital relaxation in electron excitation, where the magnitude of such effects varies 
among the different terms of a given conhguration. Furthermore, hnding numerical solutions 
that simultaneously reproduce all the the wave-function conditions can be difficult and 
iterative self-consistent treatments, e.g. Hartree-Fock, may fail to converge or yield poor 
quality results when compared with measurements (level energies and oscillator strengths). 
Calculations that approximate the wave functions by conhguration mixing tend to become 
intractable, because the collapse of the electron localizations into narrow potential wells 
gives rise to strong electron exchange interactions. Computations that employ distinct 
non-orthogonal orbitals for each conhguration or for each LS term are very difficult due to 
the large number of orbitals that need to be optimized. Additional complications arise from 
spin-orbit coupling and relativistic corrections whose ehects on calculated energy levels are 
comparable to the energy level separations. 


There is mounting observational evidence that current [Fe II] spectral models remain 
of insufficient accuracy. For example, the predicted [Fe II] line intensities in the Orion 
nebula, the archety pical H II region, disagree with observations by up to several factors (see 


Verner et ah 


2 OOOI) . In the context of extragalactic astronomy, there have been signihcant 


ehorts to use the Fe II/Mg II emission ratio as a direct Fe/Mg abundance indicator in 
quasars. According to models of cosmological nucleosynthesis, Fe enrichment trails behind 


a-element enrichment until ~ 1—2 Gyr after the initial star-formation 
many groups are activ ely trying to hnd such a point of inhexion (e.g. 


jurst; consequen 


Kurk et ah 


2007 


ly, 


Sameshima et ah 


2009|). However, most of these efforts have been inconclusive due to the 


large scatter in the Fe/Mg ratio. Uncertainties arise from the use of Fe II(UV)/Mg II as an 
abundance indicator, which are exacerbated by the fact that the classical photoionization 
models fail to account for the Fe II(A4570)/Fe II(UV) ratio by an order of magnitude; 
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therefore, Fe I 

unreliable (see 

: abundance estimates derived from using current spectral models are 

Collin-Souffrin et ah 

1980; 

Collin & JolJ 

2000; 

Baldwin et ah 

2004 

)• 


We report new calculations of Wvalues and collision strengths for the lowest 52 
even-parity levels of Fe II. The atomic data were computed using a multi-platform approach, 
where most of the state-of-the-art numerical methods of atomic physics have been used in a 
concerted effort to provide consistency checks and comparisons. Furthermore, we present 
NLTE spectral models whose predictions are benchmarked with the available astronomical 
spectra. These comparisons provide stringent tests on the quality of the atomic data. 


We also present a detailed analysis of the inherent uncertainties in the atomic data 


and their implications in N 
described by 


Bautista et ah 


.TE s 


jectral models; for this analysis, we follow the method 


is given by 


(120131) . Under steady-state balance the population of a level i 

( 2 ) 


^ki^eQki T -^ki') T 5a:*) 

Ni = - „ - ^—n = 


UeTi Qij + 1 


WT Qij + 1 


where Ug is the electron density, A^i is the Einstein spontaneous radiative decay rate 
from level k to level i, b^i is the branching ratio, Qki is the electron impact transition rate 
coefficient, and r* = ^^e level lifetime. Expressing level populations in terms 

of lifetimes and branching ratios, as opposed to using only Wvalues, has the practical 
advantage that lifetimes are generally dominated by a few strong transitions; i.e. they are 
generally more accurate than individual rates. Therefore, this level-population formalism 
gives a more clear insight into the propagation of atomic data uncertainties. 


2. Atomic structure and radiative calculations 

Realistic representations of the atomic structure are needed to obtain accurate energy 
levels, line wavelengths, and A-values, and also because such representations are the basis of 
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reliable scattering calculations. For this work we u se a com 
the pseudo-relativistic Har tree-Fock (hfr) c ode of 


Cowan 


Dirac-Fock (mcdf) code fiDvall et al 


Dination of numerical methods: 


198ll) : the Multiconhguration 


19891 ). and the scal ed Thomas-Fe rmCDirac 


central-field potential as implemented in AUTOSTRUCTURE flRadnell 


1997h . 


2.1. HER calculations 


HER uses a superposition of conhgurations approach to account for conhguration 
interactions. The code solves the Hartree-Fock equations for each electronic configuration. 
Relativistic corrections are also included in this set of equations. The radial parts of the 
multi-electron Hamiltonian can be adjusted empirically to reproduce the spectroscopic 
energy levels in a least-squares fit procedure. These semi-empirical corrections are used to 
account for the contributions from higher order correlations in the atomic state functions. 


The following conhgurations were explicitly included in the physical model; 3d®4s, 
3d®4s^, 3d^5s, 3d®4d, 3d®5d, 3(i®4p^, Sd^AsAd, 3s3p®3d’^4s, 3s3p®3d®, 


and 3s3p®3(i®4s^. ' 
HER calculation by 


his conhguration e xpansion extends the one used in the previous 


Qninet et al. 


(199a) by including 3d®4d^ and 3s3p®3d®4s^. In order 


to minimize the discrepancies between computed and experimental energy levels, the 
HER technique was used in combination with a well-known least-squares optimization 
of the radial parameters. The htting procedure was applied to S d ^As, 3d and 3d®4s^ 


with the experimental energy levels compiled by 


Sugar fc CorlissI (1198511 . In the absence 


of conhguration interaction, the 3d’^ and 3d^4s^ conhgurations are described by four 
parameters, namely the average energy F^(av), the Slater integrals F^(3d, 3d) and F^(3d, 3d), 
and the spin-orbit parameter C{3d), while for 3d®4s the exchange interaction integral 
G^(3d, 4s) is also required. In addition to these parameters, ehective interaction parameters 
such as a and /5, associated with the excitation out of the 3s and 3p subshells into the 3d, 
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are used in the fit. The average deviation between computed and experimental levels was 
found to be equal to 78 cm“h 


2.2. AUTOSTRUCTURE calculations 


AUTOSTRUCTURE fjBadnelllllQQTl. I2OI1I) computes Cl state wave functions built using 
single-electron orbitals generated from a scaled Thomas-Fermi-Dirac-Amaldi (TFDA) 
potential. The scaling factors for each orbital are optimized in a multiconhguration 
variational procedure minimizing a weighted average of LS non-relativistic term energies 
or term energies including the effects of one-body Breit-Pauli (BP) effects. Spin-orbit 
coupling and BP operators are introduced as perturbations to obtain hne-structure 
relativistic corrections. Semi-empirical corrections can also be applied to the multi-electron 
Hamiltonian, where the theoretical LS term energies are corrected in order to reproduce 
the centers of gravity of the available experimental multiplets. 


Bautistal (120081) introduced non-spherical multipole corrections to the TFDA potential 
to account for some of the electron correlation effects. This work also considered alternative 
optimization techniques of the scaling parameters for systems where the spin-orbit and 
relativistic effects are important. For the present work, we found necessary to use such 
developments but with a single variation. The modihed potential is written as 


V^(r) = H 


TFDA 


(r; A) + 


pS pSQ 

/ p{r2)slds2 + s p{s2)dS2 


( 3 ) 


where s = r/X^r, and p(r) is the electron charge density given by 


p(r) = 


27r2 1 TT 


+ 


TT^ 


+ K) - V{r) 


1 / 2 ' 


( 4 ) 


with 


Pn = 


15 2{Z-N) 


ro 


IGtt^ 


( 5 ) 

















Z being the nuclear charge and iV the e l ectron number of the system. This potential 
is different from the original in iBautistal (120081) inasmuch as being limited to a single 


dipole correction modulated by the scaling parameter A^, but more importantly, the radial 
dependence of this dipole correction is now scaled by an additional parameter \dr- 


We have performed numerous calculations with different conhguration expansions, 
starting with those in previous work then adding more conhgurations and/or using corrected 
TFDA potentials. From these computations we have arrived at the following general 
conclusions: 


(a) The inclusion of the 3s^3p^3d ^4d conhguration with a Ad pseudo-orbital is essential to 
obtain a satisfactory structure for the low even-parity levels. 

(b) The 3s^3p^3d^nl conhgurations with n > 4 have little ehect on the structure of the low 
even-parity levels. 

(c) Very large conhguration expansions tend to give inferior results than well-selected, 
concise representations. 

(d) The A-values for forbidden transitions among the low even-parity levels are very 
sensitive to the predicted energy of the 3d^As^ term. 

(e) It is difficult to reproduce the observed energy diherence between the ground state 
3d ®4s and the hrst excited term 3d ^ using the standard TFDA potential. 

(f) The magnitude of the spin-orbit correction to the energy separation between the 

3d^As and states is large; therefore, an optimization of the atomic orbitals 

based on LS energies can be misleading. 

Table [U shows a selected set of calculations in the present study. The spectroscopic 

conhgurations 3d^4s, and 3d®4s^ are common to all expansions; 3d^4d is also common 
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to all as it accounts for essential relaxation effects in configurations involving the 3d orbital. 


The hrst listed calculation is labe 
expansion of 


ed “B 


extend TFDAc”, and is built up from the original 


Bautista fc PradhanI fll996[) but taking into account additional conhgurations. 


The second calculation (“Q96+4d^”) includes the same expansion as our hfr calculation. 
In “7-con£g”, we tried to use the smallest expansion possible, comparable to the six 
configuration expansion created with the mcdf method (see next section), but with the 
addition of a Ad orbital, which could not be optimized in the MCDF method. While this 
expansion seems very small it was surprising to find that it yields some of the A-values that 
best agree with observed spectra (see Section 3). The last calculation, “BP new-TFDAc”, 
is similar to “BP extend TFDAc” but makes use of the new correlated TFDA potential 
quoted in Equation 3. 

Table [2] lists the term energies for the various calculations, showing the non-relativistic 
energies calculated in LS'-coupling as well as t he term-averaged BP e nergie s. As a reference. 


we start by looking at the Fe II expansion of 


Bautista fc PradhanI f 19961) . This model 


optimization scheme was based on the non-relativistic LS energies which are seen to be 
underestimated by 33% for the 3d'^ state and overestimated by 23% for 3d^As ^D. It 
correctly predicts the relative order of the first six LS terms, permutes the order of the 
seventh and eighth terms (dd^ and respectively) and mis-assigns the positions of 
most of the higher terms. The important 3d^As^ term is predicted to lie tenth relative 
the ground term in contrast to the observed position (12th). By including relativistic and 
spin-orbit coupling effects in this model, the predictions change considerably and for the 
worse. The term-averaged energy for the state is now overestimated by a factor 

of 2.4, and all the other terms within the 3d^ configuration also deviate farther from the 
observed energies. The relative order of the energy terms also deteriorates by including 
spin-orbit effects; for instance, the 3d^4s^ term is predicted to be as low as the eighth 
position. These large changes in the atomic structure of this model, caused by spin-orbit 
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effects, are consequences of the limitations in optimizing the Fe II atomic model on the 
non-relativistic LS energies as traditionally performed in autostructure. 


The polarized TFDA potential of 


Bautistal (120081) does not improve by itself the atomic 


model as evidenced by the predicted energies obtained with “BP extend TFDAc” (see 
Table H]). The expansion “Q96+4d ^-corr” is signihcantly smaller and is expected to be of 
inferior quality. This is conhrmed with the poor energy prediction for the state 

as well as the signihcantly higher core energies. On the other hand, relative to the ground 
term, this model predicts energies for the 3d ®4s and 3d ^ that compare with experiment as 
favorably as the other models. The last two columns in Table |2] show the results obtained 
with the new polarized TFDA potential (Equation 3) and an optimization scheme based on 
the (2j+l)-averaged energies. In order to optimize the energy of the hrst excited 3d'^ AF 
term to within 3% with respect to experiment, its non-relativistic LS energy becomes lower 
than that of the 3d ®4s ground term. This latter model is also better than the previous 
models as far as the relativistic core energy, despite the fact that the non-relativistic LS core 
energy actually looks higher. This Ending once again confirms that, in order to optimize an 
Fe II wave-function representation, the spin-orbit and relativistic energy corrections must 
be taken into account. 

Recent versions of autostructure enable the inclusion of one-body BP relativistic 
operators in the hamiltonian. Thus, orbitals can be optimized on term energies that include 
these relativistic effects, albeit missing spin-orbit splitting of fine structure levels. We 
verified that the dominant relativistic corrections to the term energies are indeed accounted 
for by the one-body operators. Hence, this new feature could have been used to optimize 


the Fe II system instead of the approach of Bautista (2008) . Though, nei t 


ler of these 


two techniques was available in the code superstructure 
previous works. 


Eissner et ah 


( I 974 J ) used in 
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Table [3] gives a complete list of the energy levels considered in this work, where the 
assigned level indexes will be the reference for the remaining of the paper. 


2.3. MCDF calculations 

Atomic descriptions were obtained within the multiconhguration Dirac-Fock (mcdf) 


framework with the GRASPO fGeneral-pur 

30se 

lelativistic Atomic Structure Package) 

code ( 

Grant & McKenzie 

1980; 

Grant et ah 

1980 

McKenzie et ah 

1980; 

Norrineton 

2004) 


Our best models included six and 12 non-relativistic configurations. The six configurations 
model included Sd^As, 3d’^, Sd^As"^, Sd^Ap, 3p‘^3d^, and 3p^3d^ (with an empty 3s orbital). 
The twelve conhguration model included 3d^As, 3d’^, 3d^As‘^, 3d^Ad, 3p^3d^, 3p^3d^As‘^, 
3p^3d^, 3s3p^3d^, 3s3p^3dHs, and 3p^3dHs. 

We found in the autostructure calculations that the inclusion of the 3d^Ad 
conhguration was important to obtain an accurate representation of the atomic structure. 
However, we were unable to obtain a fully converged Ad orbital by using the EAL (Extended 
Average Level) optimization of the 63 metastable levels. This option optimized a weighted 
trace of the Hamiltonian, the weighting factor being proportional to the statistical weights 
(2J + 1) of the levels considered. We therefore optimized all the other orbitals (Is, 2s, 2p, 
3s, 3p, 3d, As, and Ap) within an EAL optimization procedure, and obtained the Ad orbital 
in a single-conhguration mcdf calculation. 

The predicted energies of the twelve conhguration models are slightly better than from 
the smaller model, yet rather poor in comparison with experimental energies. The six 
conhguration model is much better suited for subsequent scattering calculations (see Section 
4). The average agreement between the experimental and theoretical energy levels is around 
20%. Major discrepancies are observed in levels belonging to the low-lying term 3d^As a^D 
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(around 75%). The relative ordering of the metastable even parity states is also somewhat 
problematic. When using either one of these models to compute radiative lifetimes it is 
found that the computed values disagree with all other calculations and experiments for 
most levels (up to an order of magnitude for the most sensitive levels) indicating that level 
mixing may not be properly described. 

In an attempt to improve the mcdf model, we performed a calculation using a Ad 
orbital from a multiconhguration AUTOSTRUCTURE calculation. This emploied a utility 
program to read the autostructure orbitals from a disk hie to transcribe the orbitals 
from a linear radial mesh to the GRASPO exponential prescription. However, this procedure 
neither improved the agreement between the theoretical and experimental energies nor the 
lifetimes. 


2.4. Lifetimes, branching ratios, and A-values 

As discussed in Section 1, it is convenient in NLTE spectral models to explicitly write 
level populations in terms of level lifetimes and branching ratios rather than transition 
probabilities (A-values). We discuss the radiative data from the various approximations in 
terms of these quantities, and compare them with previous results in the literature. From 
these comparisons, we arrive at recommended values and estimated uncertainties. 


Table m tabulates the present level radiative widths (Aj = t~ = A{i —>• j)) as well 
as previously published values. The last two columns indicate our recommended values and 
estimated uncertainties, which are obtain ed from the mean va l ues a nd statistical dispersion 


among all the available atomic data (see 


Bautista et ah 


2009 


20131 ). This procedure was 


applied to all levels with a few exceptions that yield grossly discrepant results with respect 
to the rest of the calculations which are then removed. This was the case of levels 21, 28, 31, 
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43, 44, 47, and 58. Another exception is the li: 


importance in the excitation of the Fe II ion flBautista et al 


etime of level 6 (3d 


7 4 


F 9 / 2 ), which is of great 


2013 ). The lifetime for this 


level is determined by a single transition, 3d'^ ^^^ 9/2 — 3d®4s ^Dg/ 2 , which is very difficult 
to render accurately due to configuration-interaction (Cl) cancelation effects. The standard 
A-value deviation in this transition is ~ 80%. Given the importance of this transition, it 
had to be treated with much more detail. 


In order to compute the rate for the highly mixed 3d^ ^^ 9/2 — ^Dq /2 transition, 

the observed transition energy must be hrst reproduced, which is only achieved by the 
model using the newly modihed TFDA potential (Equation 3). In Figured] we plot the 
calculated A{3d'^ ^P<^i 2 — 3d^4:S ^Dq/ 2 ) value vs. its predicted transition energy. These 
values are obtained from the NewTFDAc model by varying the optimization parameters in 
the potential. The A-value predicted by this model is our recommended value in Table H] 
assigning a conservative uncertainty of 30%. While we expect these models to give a 
reasonably reliable A-value for the 3d'^ '^Fg /2 — 3d^4:S ^£* 9/2 transition, it seems yield poor 
results for transitions involving higher excitation multiplets. This is apparent in Table 0] by 
the fact that the NewTFDAc model yields many outliners. 

In regards to the overall accuracy of the lifetimes, the observed dispersion between 
results of different models give an indication of how well converged the results are, thus we 
suggest that such dispersion can be use as an uncertainty indicator. In this sense, we find 
that the lifetimes for the lowest 16 levels of Fe II, which are responsible for the infrared and 
near-infrared spectra, are known within 10% or better with only a few exceptions. The 
lifetimes for higher levels, which yield the optical spectrum have uncertainties that range 
between ~10 and 30%. As to which particular atomic model is the most accurate, that is 
difficult to say from a purely theoretical point of view, thus further analysis is needed in 
view of experimental and astronomical spectroscopic information. 
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In Table [5] we present a branching-ratio sample from our various computations as 
well as published radiative data; average values and standard errors are also shown (the 
complete table is available in electronic form). It may be appreciated that the stronger 
transitions, of more practical interest, are typically those with large branching ratios, i.e. 
near unity and carry the lesser uncertainties. 

Table |5] illustrates that while branching ratios from the levels of the lowest two 
multiplets of Fe II are well known the uncertainties increase for higher levels. For example, 
in the branching ratios from level 10 there are two models that yield discrepant results. The 
“BP extend TFDAc” model yields the largest branching ratios to levels of the a “^F multiplet 
at the expense of diminishing the ratios to the ground a ®D ground multiplet. This can be 
understood from the fact that the “BP extend TFDAc” model yields a very high energy for 
the a ^F term, thus increasing the overlap of these levels with the a levels. The other 
discrepant model is the “Q96-|-4d^-corr”, which yields the smallest branching ratios to the 
a ^F levels. 

Figure [2] (left panel) shows the estimated uncertainty distribution of 387 branching 
ratios greater than 0.01. The right panel of this figure depicts the A-values uncertainty 
distribution. It may be seen that more than a quarter of all branching ratios are constrained 
to within 10% and about two thirds to within 20%. The uncertainties in the absolute 
A-values can be seen as the combined uncertainties of both branching ratios and lifetimes; 
therefore, only a small fraction (~0.02) is determined to better than 10%. A positive 
result nonetheless is that over two hfths of all transitions are in accord to within 20% and 
nearly three quarters to within 30%. We expect that the majority of observable lines and 
transitions regulating the ion population balance are sufficiently accurate. 

The complete set of A-val ues for transitions am ong the levels considered here is 


available online from AtomPy flMendoza et ah 


2014J) or by request from the authors. 
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The present data set should provide a solid platform for the modeling of Fe II spectra. 
Therefore, we argue that further work in improving the quality of the radiative data should 
concentrate on specihc transitions of observational interest rather than on the whole lot. 


3. Radiative data benchmarks 


In Section [2.41 we presented recommended lifetimes (radiative widths), radiative 
branching ratios, and d-values as a result of a critical assessment of the available data and 
from the various calculations carried out here. We also provided uncertainty estimates 
for these quantities based on the data statistical dispersion. In this section we make use 
of laboratory measurements and astronomical observations to benchmark the quality of 
the theoretical atomic data, and try to assert the reliability of the recommended data and 
their estimated uncertainties. It is important to point out here that we have used observed 
energies in the calculations of all radiative rates. 


The theoretical lifetimes presented 


in Section 12.41 are compared in Table P with 


measurements from the Ferrum Project (IHartman et ahl 


2003; 


Gurell et ah 


2009 ). It is seen 


that the lifetimes are not reproduced accurately as a whole by any of the calculations. The 
very long life of the level 43 is particularly difficult to compute accurately, and only the 
results of Deb and Hibbert (2011) agree with experiment. The results of the “new TFDA” 
model exhibit the largest differences with this experiment. This is to be expected since this 
model was optimized on transitions among the first two terms and is known not to represent 
well the highly excited levels. The scatter in the results of the various calculations owes to 
the fact that the levels are highly mixed and even small variations in the mixing coefficients 
in the different representations can lead to cancelation effects in the line strengths. Though, 
provided that all the dominant conhgurations are included in all calculations, the computed 
line strengths should all scatter around the correct answer. Thus, we suggest that by 
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averaging over the results of a number of reasonable representations of the ion one should 
arrive to reliable results. This approach seems to be well supported by Table El where our 
recommended values, resulting from averages of various results, agree with experiments to 
within the estimated error margins. 


We now compare the observed and theoretical intensity ratios between emission lines 
arising from the same upper level. In taken ratios of lines from the same upper level the 
population of that level, any dependance on the physical conditions, cancels out and the 
ratios are given by 


F{i j) _ A{i j) X{k ^ 0 _ Ki j) Kk 1) 

X “—;-^ — , ,,-X 


( 6 ) 


F(k ^ 1) A(k 1) \(i -t j) b(k ^ 1) \[i j) 

The advantage of looking at these ratios is that they depend only on A-values or, more 
specihcally, on branching ratios regardless of the physical conditions of the plasma. 
Therefore, the ratios ought to be the same in any source spectra, provided the spectra 
have been corrected for extinction. Fe II yields the richest spectrum of any astronomically 
abundant chemical species; thus, its high-resolution optical and near-IR forbidden lines are 
the best suited for the present evaluation. 137 [Fe II] lines are found in the HST/STIS 
archived spectra of the Weigelt blobs of r] Carinae. 78 [Fe II] lines are present in the 
deep ech elle spectrum (R = 30 000 ) of the Herbig-Haro object (HH 202) in the Orion 


Nebula flMesa-Delgado et ah 


2009|), and 55 [Fe II] lines have been measured in the 


X-shooter (3000 


(iGiannini et ah 


—25 0 00 A) spectrum of the jet of the pre-main-sequence star SEO-Ha 574 


201311 . 


The T] Carinae spectra are at —28° position angle (PA), and include Weigelt B and 
D positioned at 0.15" and 0.28". Six medium dispersion spectra {R = 6000—10 000) of 
the blobs have been recorded between 1998 and 2004 at various orbital phases of the 
star’s 5.5-year cycle. Two complete spectra, recorded during the broad high state and 
the several month long low state, were offset onto Weigelt D observations at PA = 68°. 
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Additional observations centered on r] Carinae recorded the spectrnm of Weigelt C offset to 
the sonthwest. Over 900 HST/STIS spectral segments were recorded of rj Carinae and the 
Homnncnlns (rednced line-by-line spectra are available onlinqj). To estimate accnrate line 
fluxes we needed to control contamination effects by stellar radiation, arbitrary continuum 
placement, and unidentihed line blends. The hrst two issues were mitigated by choosing 
slit orientations that avoided the bright central star, and by finding spectral extractions 
that minimized the continuum. To hnd blends and other extraneous features affecting the 
lines, we assumed Gaussian prohles dehned by the HST/STIS instrument response, and 
determined average centroid velocity and full-width-at-half-maximum for each species by 
fitting a sample of the stronger features with no known blends and clean continua. Then, 
all line fluxes were measured by fitting the features with Gaussians constrained by the 
parameters determined above, but allowing them to vary within the estimated uncertainties. 
We carried out five independent spectral reductions with up to four measurements of every 
observation for different spectral extractions along the GGD and different assumptions 
about the continuum and noise levels. Our visible spectrum measurements showed that 
line centroids and widths can be constrained to within 2 kms“^; therefore, blends and 
contamination that affected line peak by more than 0.03 A could be readily identified. Line 
fluxes of strong, unblended features were measured with an accuracy of better than 10%. 
For other less certain lines, we at least obtained robust error bars. 


Measured line fluxes in rj Garinae must be corrected for extinction, and this affects 
the comparison between observed and theoretical line ratios. While the extinction curve 
toward the Weigelt blobs and other regions of the nebula are not well understood, most 


spectroscopic evide nce suggests t 


and Ry = 3.1 (see 


lat th e extinction curve is well described by Ay = 2.1 


Bautista et al 


2011 


, for a discussion). Here we explore the extinction 


^http://etacar.umn.edu/ and http://archive.stsci.edu/prepds/etacar/ 
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magnitude and selective extinction, and find that the above parameters yield the best 
agreement with the theoretical line ratios. Moreover, it is found that even large variations 
in the adopted visual extinction magnitude would not change the our conclusions regarding 
the Fe II A-values choice. 


Comparisons of measured line ratios show significant systematic errors that are easily 
overlooked, and by analyzing multiple measurements of the same line ratio, we attempt to 
minimize such systematic errors and gain insight on the observational accuracy. Regarding 
the Tj Carinae spectra, we find 107 reasonably well-measured line ratios that are defined as 


ratio = 


max(Fl, F2) 
min(Fl,F2) ’ 


(7) 


where FI and F2 are the measured fluxes of two lines from the same upper level. The 
minimum of the two fluxes is in the denominator such that the line ratios are unconstrained, 
and equally weighted when compared with the theoretical expectations. Figure [3] shows 
a line-ratio sample from various observed spectra as well as from multiple theoretical 
determinations, where some of the most common traits are illustrated. A regular feature is 
the scatter in the observed values which greatly exceeds the estimated individual errors. 
Also the scatter in the measured line ratios often exceeds that of the theoretical predictions. 
In some cases (see the bottom-left panel of Figure [3]), there are systematic differences 
between the measured ratios in rj Carinae, SEO-Ha 574, and HH 202. Note that we only 
measured the lines in rj Carinae, while for SEO-Ha 574 and HH 202 we adopt the published 
line intensities. Thus, the fact that there are discrepancies between the ratios from the 
latter two sources is further evidence of unaccounted systematic errors in the measurements, 
independent of our own work. In five cases, the differences between observation and theory 
are in excess of the scatter of the different methods (see the bottom-right panel of Figure E]). 
(Plots of all 106 line ratios extracted from observations are available from the authors upon 
request.) The results of these comparisons should serve to warn researchers against the 
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temptation to try to derive atomic parameters from a single set of observed spectra. 


The adopted line ratios for the present study are the mean observed values and their 
uncertainties are given in terms of their standard deviations. These line ratios are then 
compared with theory and our recommended branching ratios. It is worth mentioning that, 
while the ratios determined from ^4-values and branching ratios are equivalent, the error 
margins determined from the latter are smaller and thus preferable. Table [7| shows reduced 
statistical indexes for this comparison which are dehned by the expression 


— 2 / ^ (-Rob ~ -Rth)^ 

^ =\In2^' 


( 8 ) 


N ^ (SRl^ + SRi^) • 

where Roh and i?th are the observed and theoretical ratios, respectively, and SR is the ratio 
uncertainty. The reduced quantihes the agreement between theory and observations: 
perfect agreement (within the adopted uncertainties) gives a value of 1.0; discrepancies 
beyond the stated uncertainties lead to > 1; and a reduced < 1 indicates that the 
adopted uncertainties have been overestimated. 

The hrst thing to note from Table [7] is that the observational uncertainties by 
themselves are in general insufficient to account for the discrepancies between theore tical 


line r atios and observations. We see that the recent CIV3(DH11) flDeb fc Hibbert 


201lh calculation is marginally better than the previous computation by 


Quinet et ah 


fjl996l) . The present NewTFDAc results look someone better than previous calculations. 
This seems contrary to the lifetimes comparisons of Table El This is because absolute 
A-values and radiative line widths are a lot more difficult to compute accurately. 
Thus, it is a reasonable practice to measure lifetimes experimentally and use them 
to correct A-values, while keeping the theoretical branching ratios. The ‘7-con£g’ 
calculation, which employs the smallest expansion considered here, does yield the best 
agreement with observed spectra. Our recommended values seem to be a significant 
improvement over the older theoretical data, and once the theoretical uncertainties 
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are taken into account, reasonable agreement with observations is reached. It may be 
appreciated that the branching ratios uncertainties are smaller than the uncertainties 
in the radiative rates, and if the latter are adopted = 0.87, thus indicating 
that the uncertainties are overestimated. The larger departures between theory and 
observations are associated to hve line ratios: A4458.5/A4515.5 (1(33—2)//(33—3)); 
A4320.1/A4147.6 (J(39-8)/J(39-6)); A4373.0/A4147.6 (J(39-9)/J(39-6)); A16769/A15335 
(/(ll—7)//(l— 6 )); and A4515.5/A4775.3 (1(33—3)//(33— 6 )). Not surprisingly, four of these 
problematic ratios involve level 6 (a ^^ 9 / 2 ) which is very difficult to represent theoretically. 
If these hve ratios are removed, the reduced = 1-07. It is unclear if the remaining 
line-ratio discrepancies arise from observations or theory; nonetheless, as shown in Figure HI 
the overall agreement between theory and observations is satisfactory. 


4. Electron Impact Collision Strengths 

For the different Fe II target expansions described above, collision strengths have been 
computed by two different methods, specihcally R-matrix-|-ICFT (RM-I-ICFT) and DARC. 
In the case of the small 7-con£guration expansion it was also possible to do a Breit-Pauli 
calculation (BPRM), which compared very well with the RM-I-ICFT results. We have also 
carried out calculations with different close-coupling expansions in the stg2 R-matrix 
code, with and without energy-corrected excitation thresholds at the Hamiltonian matrix 
diagonalization in stg3. This allows us to estimate the sensitivity of the collisional data 
to target representation as well as to the details of the scattering calculation. Our runs 
include 20 continuum orbitals for each angular momentum in the close-coupling expansion, 
and partial waves up to L = 10 in LS'-coupling and J = 14.5 in JJ-coupling. Collision 
strengths are sampled at 10 000 equally spaced energy points up to the highest excitation 
threshold, with a much coarser mesh beyond it (3xthreshold). 
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For our DARC calculations the target orbitals and energy lev els were generated 


using the Dirac-Ha rtree-Fock atomic structure program GRASPO flDvall et al 


1989 


Parpia et ah 


19961) . We employed only a minimal set of 6 configurations, including 


3(f, 3p^3d^, 3p^3d^, 3d®4s, 3(i®4p and 3d^4s^ for a total of 329 levels. The first 20 levels were 


parallel Dirac R-mat 

of the original codes 

rix programs {. 

Ballance et al. 

2006) 

, which merges of modified versions 

developed bv ^ 

Vorrington 

(200 

4) with our suite of parallel Breit-Pauli 

R-matrix programs ( 

Badnell et al. 

2004; 

Mitnik et al. 

2003 

). The size of the R-matrix 


’box’ was 13.3 an and we employed 12 basis orbitals for each continuum-electron angular 
momentum. This was more than sufficient to span electron energies up to 6.0 Ryd (81 eV). 
As the model considers only low temperature modeling, we only included partial waves up 


to J= 10 with a top-up procedure to account for higher partial wave contributions 


Burgess 


(1197411 . which were minimal. The cross sections spanned the energy range from the ground 


state to just over 1 Ryd, with an energy resolution of 10,000 points. 


We computed thermally averaged effective collision strengths 

poo 

Tij= / Vtij{Ej)exp{-Ej/kTe)d{Ej/kTe), (9) 

Jo 

where Qij{Ej) and Ej are the collision strength and incident electron energy relative to the 
jth level, respectively, Te is the electron temperature, and k the Boltzmann constant. 

Effective collision strengths at 10^ K from different calculations are compared in 
Table [HI Only a sub-set of transitions from the ground level are shown, but the complete 
table for transitions among all the 52 levels is available electronically. Also, we do not 
show the results of all different calculations, but only the most representative ones. The 
mean and standard deviation of the collision strength for each transition are also listed. 
The most relevant comparisons involve excitations from levels of the ground and hrst 
excited multiplets since they dominate the whole spectrum, while most other collisional 
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transitions among excited levels have very little impact on the resnltant spectrum. Also, the 
comparison of effective collision strengths at 10^ K is perhaps the best way to make sense of 
the collision strength variations since this is the temperature where Fe II is most frequently 
found. In comparisons at higher temperatures, the differences tend to be smaller because 
the contribution of near-threshold resonances to the collision strengths is less conspicuous. 
At lower temperatures the differences are enhanced because the near-threshold resonance 
structures dominate, but their exact position and interference patterns are uncertain. 

In fact, it is very difficult to provide effective collision strengths at temperatures below 
~ 5000 K with reasonable conhdence. 


The hrst four columns of collision strengths in Table [8] sh ow the results of RM -I-ICFT 


calculations using a target with the same conhgurations as in 


Quinet et ah 


dlOOd). We 


refer to these calculations as ’Q-I-RM’. For this target, we look at the results when the 
theoretical target energies are used through the calculation (Q-|-RM-ns), the level energies 
are shifted in the H matrix to experimental values (Q-RM-shift), the radius of the R-box 
is limited to 8 a.u. (Q-|-RM-RA=8), and the radius of the R-box is extended to 14.5 a.u. 
(Q-|-RM-RA=14.5). The next set of results corresponds to an RM-I-ICFT calculation using 
the ’7-con£g’ target described in Section 2. For this calculation we chose not to correct the 
target energies in stg3. The sixth set of effective collision strengths shown is from our best 
DARC calculation. 


In addition to the calculations described above we did other calculations employing 
larger target expansions, but the results were neither signihcantly different from those 
published before nor yielded predicted spectra that compare favorably with experiment 
(see the next section). In comparing the results of various calculations we find significant 
systematic differences in the collision strengths for transitions among levels of the 3d ®4s 
ground multiplet. For instance, for the excitation from the ground level to the hrst 
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excited level, most of our current calculations give T(lO^i^) of about 2, while the darc 
calculation and previously published works give T(lO^i^) of about 5. The differences 
are found not to depend on the conhguration expansion. Instead, the effective collision 
strengths are greatly enhanced when the excitation threshold are shifted to experimental 
values. This is illustrated in Figure |5l where we plot the DARC collision strengths with 
and without energy corrections. It can be seen that by shifting the target energies the 
Rydberg series of resonances seem to be pushed against the excitation threshold, then 


the calculations of 


Zhane; and Pradhan 

(1995 

) and 

Ramsbottom et ah 

(2007) 


we also show the collision strengths using the “NewTDFAc”, which yields very accurate 
energies for the lowest Fe II terms without corrections. It is found that the “NewTDFAc” 
target does not give the large packs of resonances found in other calculations when the 
target energies are shifted to experimental values. While shifting threshold energies to 
experimental values is a technique generally used to improve the cross sections, such shifts 
are expected to be small. However, most target representations of Fe II yield very poor 
energies for the hrst few excited terms of the system, thus large corrections are needed. 
While it is unclear whether the ‘NewTDFAc’ target yields accurate collision strength for 
these transition, it does show that previously published collision strengths for the ground 
multiplet are much more uncertain than previously thought. 


By comparing the Q+RM-ns and Q+RM-shift results that shifting the energies to 
experimental values has a very large effect for several transitions to highly excited levels. 
There are two reasons for this. Firstly, shifting the target energies also shifts the Rydberg 
series of resonances, particularly those resonances that lie near excitation threshold. 
Secondly, excited metastable levels with equal j quantum numbers can be strongly mixed, 
and such mixing is greatly affected their energy separation. Moreover, when the energy 
shifts result in changing the relative order of such levels this can lead to large errors in the 
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couplings that yield the resonances on the cross sections for such levels. 

By comparing columns two, four, and hve one hnds that changing the size of the 
R-matrix box has a smaller effect on the cross sections than shifting the target energies, 
yet the effect is still sizable for a few transitions. These effects should also be considered in 
assessing the uncertainty in the cross sections. 

It should be stressed that all calculations generally agree on background cross sections, 
and the differences are mostly due to the resonance structures. Furthermore, we take into 
account all results, and their statistical dispersion is assumed an accuracy indicator. 

Regarding the overall error distribution, it is found that, for excitations from the 
lower nine levels of the ion that dominate the entire spectrum, uncertainties in the range 
of 10-20% are the most frequent. This is closely followed by those in the 20-30% range. 
Moreover, most of the effective collision strengths have uncertainties less than 50%, with 
a tendency for the smaller collision strengths to be poorer than the weaker ones. If the 
complete collisional inventory is considered, the effective-collision-strength discrepancies 
vary widely, reaching factors of 2 or more in some cases. Despite these bulk uncertainties 
in the collision strengths the strong transitions that dominate the spectrum under typical 
nebular conditions seem to be reasonably well known, as will be discussed in the next 
section. 

It should be noted that because the largest atomic data uncertainties are in the collision 
strengths it is expected that modeled spectra will be more accurate as the electron density 
increases, while larger error are expected in models of very diluted plasmas (see Bautista et 
al. 2013). 

The current data set establishes a hrm foundation from which specihc transitions of 
interest can be studied individually in order to reduce their current uncertainties. However, 
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it should not be assumed de facto that the larger expansions yield the more accurate results 
unless they are proven to be fully converged. 


5. Comparison with observed spectra 


In order to assess the collisional data quality, we have constructed collisional excitation 
models with the recommended d-values from Section 12.41 and the effective collision strength 
data sets introduced in Section 01 The predicted emission spectra for each of these models, 
for several electron temperat ures and densities, are com pared with the spectra of the 


Herbig-Haro object HH 202 


SEO-Ha 574 f Ciannini et ah 


Mesa-Delgado et ah 


20091) and the pre-main-sequence star 


20131 ). All [Fe II] lines measured in these objects are intended 


in the comparisons, i.e. 78 and 55 respectively. We make no attempt to compare with the 
[Fe II] spectrum of the Weiglet blobs of rj Carinae because the strong UV Fe II emission 
from this object indicates that there is signiheant fluorescent exci tation. Atomic rate e rror 
propagation in predicted line emissivities has been examined by 


Bautista et ah 


(1201311 : 


for the present comparisons, we adopt the statistical scatter of the A-values and effective 
collision strengths at 10^ K as estimates of the atomic data error for each transition. Instead 
of comparing absolute line intensities which depend on ion abundances, we normalize the 
observed lines to the sum of all line fluxes and the theoretical emissivities to the sum 
of the line emissivities of all lines in the observed spectrum, and we vary the electron 
temperature and density in each model to get the best possible accord with observations. 
The theory-experiment correlation is expressed in terms of the reduced 

Table [9] shows the model physical conditions and reduced values that best match 
the observed HH 202 spectrum, where each model implements a specihe set of effective 
collision strengths. As gauged by the reduce d index, it may be seen t hat the model with 


the effective collision strengths compiled by 


Bautista &: PradhanI (119981) leads to a fairly 
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poor correlation with the observed spectrum. 


Model RH07, that employe es the effective collision strengths of 


Ramsbottom et al. 


(1200711 and 


Ramsbottom 


(120091 1 fairs even worse than previous models when compared 


with optical spectra. 


The collision strengths computed with the fully relativistic DARC code yield the worst 
agreement with observation. This is to be expected, since it was not possible to construct a 
self-consistent 4d orbital for this target, that accounted for important relaxation effects of 
the spectroscopic 3d orbital. 


The RM-I-ICFT calculations shown in Table 9 use the same configurations as the 
‘New-HFR’ for the radiative calculations. This target expansion is used for two calculations 
with 63 and 114 levels in the close coupling expansion. The two results yield line intensities 
that agree significantly better with observations than previous calculations. 


The observed spectrum is best reproduced with the collision strengths of the ‘7-config.’. 
With respect to the HH 202 spectrum, this collisional data set yields a reduced of 1.01 
for Te = 11500 K and Ug = 6.56 x 10^ cm“^; i.e. a nearly perfect match. It is interesting 
that the Fe II spectrum points to a n elec tron density that is roughly four times the density 


diagnosed by 


Mesa-Delgado et al. 


(120091) with higher ionization species. This is consistent 


with an observed shocked region where the Fe II emission arises from the shock front itself. 
Figure [6] depicts a comparison between predicted and observed spectra. 


As good as the agreement with the optical spectrum of HH 202 is, it does not constrain 
the collision strengths among the lowest 16 levels of Fe II, which yield infrared lines. Thus, 
we look at the spectrum of SEO-Ha 574, which combines optical and near-IR lines. The 
between 1 to 2 m, includes 14 lines from de-excitations from the 3d®4s a multiplet to 
the 3d®4s a ®D and 3d^ a ^F terms. We find that the best agreement with observation 
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is found when combining the RH07 collision strengths for the lowest 16 levels and the 
‘7-con£g’ collision strengths for the higher levels. This model yields a reduced of 1.32 
fort Te = 10 000 K and Ug = 4.4 x 10^ cm“^. It is apparently clear that the RH07 model 
yields more accurate collision strenghts than our preferred model for transitiona among the 
hrst three terms of Fe II, and particularly between the a and a terms of the 3d®4s 
configuration. This may be due to the way the target wavefunctions were optimized in 
RH07. Yet, this same accuracy does not seem to be maintained through higher levels. 

The agreement with SEO-Ha 574 is worse than for HH 202 due in part, we believe, to 
the scantier number of spectral lines (only 53 lines) and to underestimated observational 
uncertainties. In published spectra, it is claimed that nearly a third of the lines are 
accurate to better than 5% and almost two-thirds to better than 10%. Given the lower 
resolution of the SEO-Ha 574 spectrum, relative to that of HH 202, the trial uncertainties 
in the measured fluxes could be roughly twice as large as those claimed, in which case the 
agreement with the theoretical prediction would be comparable to that found for HH 202. 


5.1. Emission spectra: IR and optical 

We examine all lines longwarth of 3 /im since different instruments cover different 
segments of the mid-IR and IR regions, the most important [Fe II] lines in increasing 
wavelength order being at: 5.34 /im {Sd^As ^Dg /2 — 3d'^ "^^ 9 / 2 ; 1 ~ 6); 17.93 /im 
(3d^ ^Fq /2 - 3d^ 6 - 7); 22.92 fxm {3dHs - 3dHs ^Fg/a; 10 - 11); 24.51 pm 

(3d^ ^F 7/2 - ^Fg/a; 7 - 8); 25.98 pm (3d®4s ^Fg/a - 3d®4s ^Fr/a; 1 - 2); 35.34 pm 
(3d®4s 6F7/2 - 3dHs ^Fg/a; 2 - 3); and 35.77 pm {3dHs ®Fg/2 - 3d^As ^Fg/a; 8 - 9). 

Figure [7] shows emissivity line ratios among these lines as a function of electron density 
for temperatures between 5000—15 000 K. These ratios can be used as diagnostics of 


28 


electron density in the range of 10^—10® cm“^. Ratios that involve the 5.34 /rm line carry a 
signihcant error, particnlarly toward the high densities dne to the 30% nncertainty in the 
lifetime of the 3d ^ ‘^^9/2 level. All other line ratios are well constrained and shonld provide 
reliable diagnostics. The j(35.34 /im)/j(25.98 /rm) ratio is particnlarly relevant becanse it 
is essentially invariant with temperatnre; therefore, this ratio in combination with any of 
the other IR ratios would constrain density and temperature in the [Fe II] emitting region. 

There is a good number of [Fe II] lines in the near-IR region (1—3 /im) that 
originate from radiative transitions involving levels within the multiplet. 

Having more than one line from the same upper level is useful as they can be used as a 
dust-extinction diagnostic. The stronger lines in this part of the spectrum are at 1.257 fim 
{3d^4:S ^Dg /2 — 3d^4:S 1 - 10) and 1.644 /xm {3d'^ ^^ 9/2 — Sd^As 6 - 10). The 

intrinsic j(1.644 /im)/j(1.257 /xm) ratio is estimated to be 0.80 with an accuracy of ~ 20%, 
but extinction due to dust would increase this ratio in observed spectra. 

The near-IR line ratios are insensitive to temperature variations at 10^ K, but they 
can be used as density diagnostics in the range 10^—10^'^ cm“^. Such line ratios reach the 
high-density limits beyond 10 ® cm“®, but near those limits, the ratios are of little use due 
to the current uncertainties in the radiative branching ratios as illustrated in Figure [HI 
Moreover, observations of dense nebulae in the near-IR could be very useful to constrain 
the high-density limits of various diagnostics and, hence, the transition branching ratios. 

The stronger lines in the red part of the spectrum (7400—10 000 A) result from 
de-excitations of the 3d ^ and levels. There are four particularly strong lines whose 
strengths are accurately determined both theoretically and observationally, and are in good 
agreement in the spectrum of HH 202 : 8616.8 A {3d'^ '^Fq /2 — 3d’^ ^^ 5 / 2 ; 6 - 14); 8891.8 A 
(3d^ ^Fy 2 - ^P 3 / 2 ; 7 - 15); 7155.2 A (3d^ ^F ^/2 - 3d^ ‘^G^i 2 ] 6 - 17); and 7452.6 A 

(3(i^ ^Fi /2 — 3d'^ ^ 1 ^ 7 / 2 ; 7 - 17). From these lines four diagnostic ratios can be determined 
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that cover the range 10 ^— 10 ^ cm for densities higher than 10 ^ cm the ratios are 
cnrrently too nncertain for any diagnostic pnrposes. Usefnl ratios are shown in FignreO 

The visible [Fe II] spectrnm (4000 < A < 7400 A) is very rich, exhibiting lines from 
excited levels ~ 2.5—3.8 eV above the gronnd level. We have selected the seven strongest 
and most reliable visible lines for diagnostic ratios: 5527.4 A ^^ 7/2 — 3d’^ ^-^ 5 / 2 ; 7 - 27); 
5158.8 A (3d^ ^F 7/2 - 3d 64 s ^idi 3 / 2 ; 7 - 26); 5261.6 A (3d^ ^^ 7/2 - 3d®4s ^idii/ 2 ; 7 - 29); 
5111.6 A (3d^ ^^ 9/2 - 3d®4s ^idii/ 2 ; 6 - 29); 5333.6 A (3d^ - 3d®4s ^idg/a; 8 - 30); 

5220.0 A (3d^ ^^ 7 / 2 -3d®4s ^idg/a; 7-30); and 4745.5 A (3d ^ ^Fg/a - 3d ^4s ^Fg/a; 6-34). 
The most usefnl line-ratio diagnostics are shown in Figure [TUI which can be used in a wide 
density range ( 10 — 10 ® cm“^) and are mostly insensitive to temperature variations around 
~ 10 000 K. 

The best temperature diagnostics are obtained from line ratios of optical, red, or 
near-IR lines. This is illustrated in Figure ITT] using the stronger lines in each part of the 
spectrum. These ratios can be very useful in the analysis of spectra from instruments such 
as the X-Shooter spectrograph that covers visible and near-IR ranges simultaneously. 


5.2. Absorption spectra: UV transitions 


UV absorption spectra of bright sources, e.g. AGN and GRB, often exhibit Fe II 
troughs. In some objects it is possible to measure troughs from excited levels together with 
the resonant transitions, and this offers the possibility to use measured 


column densities as density and temperature diagnostics (see, for example. 


Dunn et ah 


2010 


, and references therein). At typical temperatures of ~ 10^ K, the most populated 


levels and those that have been observed in absorption arise from the 3d®4s ®F multiplet 
(with energies 0.0, 384.8, 667.7, 862.6, and 977.0 cm“^) and the high-multiplicity levels 
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of the and 3(i®4s multiplets (at 1872.6, 2430.1, and 7955.3 cm”^). Fignre [12] 

shows column density ratios for these levels as a function of electron density for three 
different temperatures. It can be seen that these column density ratios are distinct density 
diagnostics for densities up to 10® cm 


It is noted that the current results resolve the discrepancy found by 


Dunn et ah 


fj2010l) in the column density from the 1873 cm ^ level (a '^'Pq/ 2 ) in the spectrum of the 


Fe 

of 


jOwBAL quasar SPSS J03 


8-0600. In that work, it was found that the collision strengths 


Bautista fc Pradhanl (119981 ) reproduced the observed column density for this level for 


a density ~ 0.5 dex lower than the estimates from other Fe II and Si II levels. Some 
fluorescence effects were put forward as possible explanations for this effect; however, the 
current atomic data, on there own merits, yield column densities much more consistent with 
observations than previous models. 


6. Discussion and Conclusions 

We present a complete spectral model for the Fe II ion comprising the lowest 52 
metastable levels of the ion. It accounts for essentially all dipole forbidden lines of the ion 
in the optical and IR spectral regions, as well as the column densities of absorption troughs 
in the UV. This model is the result from extensive revisions and calculations of atomic 
parameters, namely dipole forbidden A-values and electron impact collision strengths. 

For these calculations we employ several different state-of-the-art theoretical methods, 
and compare the results with previously published data. These comparisons allow us to 
estimate the uncertainty in each atomic rate, and propagate them through spectral model 
predictions. 

A general conclusion that can be derived from all the different calculations is that—for 
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such a complex atomic system as [Fe II] and the very large number of energy levels and 
transitions involved in the spectrum—not a single calculation can achieve convergence so as 
to provide all accurate atomic parameters at once. Moreover, very large atomic expansions 
tend to give worse overall results than small, well-controlled, and optimized expansions. 

Furthermore, we employ a number of astronomical observations with exceptionally rich 
[Fe II] spectra as benchmarks for the atomic data. We hnd very good agreement between 
observations and our recommended data within the estimated uncertainties derived from 
measurements and calculated values. Moreover, our spectral model is able to predict optical 
emission spectrum of over 100 lines of the HH 202 object in the Orion nebula in nearly 
perfect statistical agreement with observations. 

Our atomic model is then used to explore the most important [Fe II] line diagnostic 
ratios in the optical, near-IR, and IR regions. We also present useful column density 
diagnostic ratios in the UV. 

The atomic data, with estimated uncertainties, from this work is available in Table 4 
(radiative life times). Table 10 (in electronic form; branching ratios). Table II (in electronic 
form; ratios from lines from the same upper level in rj Carinea), Table II (in electronic 
form, effective collision strengths), and Table 12 (in electronic form, estimated effective 
collision strength uncertainties). 
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Fig. 1.— Theoretical A-values as a function of energy separation for the 3d'^ — 

3d^4s ^Dgi 2 transition calculated with the modihed TFDA potential model (Equation 3). 
The horizontal dotted line depicts the experimentally observed energy separation. 
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Fig. 2.— Branching-ratio and A-values error distributions for dipole-forbidden transitions 
among metastable levels of Fe II. 
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Fig. 3.— Emission-line ratios among transitions from the same upper level. The first nine 
points from left to right (’77 Carinae’) are the results from our measured intensities in the 
HST/STIS spectra of the Weigelt blobs of 77 Carinae. The tenth point (’TT’) is the measured 
ratio in the X-shooter spectrum of SEO-Ha 574. The last point to the right (’<R>’) depicts 
the average of all measurements and uncertainties given by the standard deviation. The 
horizontal lines represent the predictions from several different computations of A-values. 
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Fig. 4.— Comparison of theoretical and observed intensity ratios for lines arising from the 


same upper level. 
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Fig. 5.— Collision strengths for the 3d®4s ®D 9/2 (ground level) to the ®D 7/2 (hrst excited 
level). The top panel shows the collision strength from a LS calculation with algebraic hue 
structure splitting using the ‘NewTDFAc’ target, which gives accurate threshold energies. 
The second and third panels show the collision strengths from our DARC calculations, keep¬ 
ing the theoretical energies and shifting the thresholds to laboratory energies, respectively. 
The positions of the 3d^ ^F and 3d®4s thresholds are indicated in each panel by vertical 
dashed lines 
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Fig. 7.— IR density diagnostic line ratios. Line-ratio errors resulting from the atomic data 
are indicated by the dotted lines. The ratios are calculated at temperatures of 5000 K 
(lowest, red lines), 10 000 K (middle, black lines), and 15 000 K (upper, blue lines). The 
different ratios are: R1 = j(35.77/rm)/j(25.98/im); R2 = j(35.34/rm)/j(25.98/im); i?3 = 
j(24.52/im)/j(25.98/im); RA = j(22.89pm)/j(25.98pm); R5 = j(17.93pm)/j(5.34pm); 
and i?6 = j(17.93pm)/j(25.98pm). 
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Fig. 8.— Near-IR density diagnostic line ratios at 10 000 K. Line-ratio errors resulting 
from the atomic data are indicated by the dotted lines. The different ratios are: R1 = 
j(1.600/im)/j(1.644/im); R2 = j(1.534/im)/j(1.644/im); i?3 = j(1.294/im)/j(1.257/im); 
and i?4 = j(1.279/im)/j(1.257/rm). 
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Fig. 9. — Red density diagnostic line ratios at 10 000 K. Line-ratio errors resulting from 
the atomic data are indicated by the dotted lines. The different ratios are: R1 = 
j(7155.2 A)/j(8616.8 A); R2 = j(7452.6 A)/j(8616.8 A); R3 = j(7155.2 A)/j(8891.8 A); 
and RA = j(7452.6 A)/j(8891.8 A). 
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Fig. 10. — Optical density diagnostic line ratios at 10 000 K. Line-ratio errors re¬ 
sulting from the atomic data are indicated by the dotted lines. The different ra¬ 
tios are: R1 = j(5261.6 A) /j(5333.6 A); R2 = j(5158.0 A)/j(5333.6 A); R3 = 
j(5527.4 A)/j(5333.6 A); R4 = j(5158.0 A)/j(5261.6 A); R5 = j(5527.4 A)/j(5261.6 A); 
and R6 = j(5158.0 A)/j(5527.4 A). 
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Fig. 11. — Temperature diagnostic line ratios combining optical and red or near-IR lines. 
The ratios are plotted for densities of 100 cm“^ (black lines), 10^ cm“^ (green lines), 
and 10® cm“® (blue lines). Line-ratio errors resulting from the atomic data are indi¬ 
cated by the dotted lines. The different ratios are: R1 = j(5261.6 A)/j(8616.0 A); R2 = 
j(5158.0 A)/j(8616.0 A); R3 = j(5261.6 A)/j(7155.2 A); R4 = j(5158.0 A)/j(7155.2 A); 
R5 = j(5261.6 A)/j(1.2567 /xm); and R6 = j(5158.0 A)/j(1.2567 pm) 
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Fig. 12.— Excited level column density (relative to that of the ground level) as a function 
of electron density. Results for temperatures of 5000 K, 10 000 K, and 15 000 K are listed 
from left to right. The top panels show the results for the four excited levels of the 
multiplet, starting with the level with E = 384.8 cm“^ (top curve) and ending with that 
with E = 977.0 cm“^ (lowest curve). In the lower panels we show the relative column 
densities of levels with E = 1872.6 cm“^ (top curve), E = 2430.1 cm“^ (middle curve), and 
E = 7955.3 cm“^ (lowest curve). 





















Table 1. autostructure configuration expansions for Fe II 


Label 

Configuration Expansion 

Spectroscopic 

3s'^3p^3d'^, 3s^3p^3d^4:S^ 

BP extend TFDAc 

3s^3p^3d^Ad, 3s^3p^3d^Ap, 3s^3p^3d^AsAp, 3s^3p^3d'^As, 3s'^3p^3d^, 3s^3p^3d^, 
3s3p^3d'^As, 3s3p®3d’^4p, 3s3p®3d®, 3p®3(i®4s, 3p®3d’^4s4p, 3p®3d®, 3p®3d’^4s2 

Q96+4d2 

3s^3p'^3d'^Ad, 3s3p®3(i’^4s, 3s3p®3(i®, 3s'^3p'^3d^Ap'^, 3s^3p'^3d^Ad^, 3s^3p^3d^AsAd, 
3s^3p<^3d<^5s, 3s‘^3p*^3d^M 

7-config 

3s^3p®3d^454d, 3s3p®3d^4s, 3s3p®3d® 

BP new-TFDAc 

3s23p®3(i®4d, 3s'^3p^3d^AsAp, 3s'^3p^3d^Ap, 3s'^3p^3d'^As, 3s'^3p^3d^, 3s^3p^3d^, 

3s3p®3d'^4s, 3s3p®3(i’^4p, 3s3p®3d®, 3p®3d®4s, 3p^3d'^AsAp, 3p®3d®, 3p®3d'^4s2 


Note. — Spectroscopic configurations give rise to the levels of interest and are common to all expansions. 











Table 2. Comparison of autostructure energies with experiment 


Term 

Expt 

BP 

BP exted TFDAc 

Q96+4d 

^-corr 

NewTFDAc 

LS 

(JJ) 

LS 

(JJ) 

LS 

(JJ) 

LS 

(JJ) 

3d®4s ®D 

0.000 

0.000[01] 

0.000[01] 

0.000[01] 

0.000[01] 

0.000[01] 

0.000[01] 

0.001 [02] 

0.000[01] 


0.018 

0.013[02] 

0.034[02] 

0.062[02] 

0.077[02] 

0.034[02] 

0.058[02] 

0.000[01] 

0.019[02] 

3d64s *D 

0.072 

0.089[03] 

0.093[03] 

0.091[03] 

0.090[03] 

0.087[03] 

0.093[03] 

0.085[03] 

0.084[03] 

3d'^ 

0.120 

0.124[04] 

0.148[04] 

0.164[04] 

0.177[04] 

0.158[04] 

0.180[04] 

0.110[04] 

0.127[04] 

3d'^ 

0.143 

0.154[05] 

0.179[05] 

0.199[05] 

0.214[07] 

0.175[05] 

0.199 [05] 

0.141[05] 

0.160[05] 


0.165 

0.167[06] 

0.192[06] 

0.212[07] 

0.226[08] 

0.208[06] 

0.231 [07] 

0.155 [06] 

0.173[06] 

3d'^ 

0.184 

0.208[08] 

0.234[10] 

0.247[12] 

0.271[14] 

0.220[08] 

0.244[10] 

0.194[08] 

0.214[08] 

3d'^ 

0.186 

0.187[07] 

0.214[07] 

0.232[09] 

0.249[11] 

0.230[10] 

0.255[11] 

0.174[07] 

0.195 [07] 

3d®4s 4p 

0.191 

0.229[11] 

0.235[11] 

0.224[08] 

0.226[09] 

0.224[09] 

0.232[08] 

0.230[10] 

0.230[10] 

3d®4s 

0.192 

0.208[09] 

0.228[09] 

0.211[06] 

0.212[06] 

0.212[07] 

0.218[06] 

0.224[09] 

0.224[09] 

3d®4s ■‘P 

0.204 

0.291[15] 

0.249[12] 

0.232[10] 

0.191[05] 

0.234[11] 

0.239[09] 

0.245[12] 

0.245[12] 

3d®4s2 65 

0.209 

0.214[10] 

0.218[08] 

0.243[11] 

0.232[10] 

0.425[16] 

0.403[16] 

0.236[11] 

0.234[11] 

3d®4s 

0.231 

0.275[12] 

0.280[13] 

0.256[13] 

0.258[12] 

0.260[12] 

0.267[12] 

0.276[14] 

0.278[14] 

3d®4s 2p 

0.235 

0.281[14] 

0.287[15] 

0.277[15] 

0.279[15] 

0.275[14] 

0.284[14] 

0.280[15] 

0.280[15] 

3d®4s 2R 

0.235 

0.276[13] 

0.283[14] 

0.267[14] 

0.266[13] 

0.264[13] 

0.273[13] 

0.274[13] 

0.273[13] 

3d®4s 2p 

0.246 

0.291[16] 

0.298[16] 

0.283[16] 

0.284[16] 

0.283[15] 

0.290[15] 

0.286[16] 

0.290[16] 

Core Energy 


-2324.32 

-2541.73 

-2524.32 

-2541.70 

-2523.77 

-2540.17 

-2523.14 

-2542.70 


Note. — Energies are given in Ryd. Numbers in square brackets indicate the relative term positions. 







Table 3. Energy levels of the 3(i®4s, 3d^, and 3(i^4s^ configurations 


Index 

Configuration 

Level 

Energy* (Ry) 

Index 

Configuration 

Level 

Energy* (Ry) 

1 

3(l®(®D)4s 


0.0000000 

27 

3d®(®P2)43 

^Pl/2 

0.2042136 

2 

3(l®(®D)4s 

®^7/2 

0.0035065 

28 

3d® (377)45 

'^Pl3/2 

0.1936589 

3 

3(l®(®D)4s 


0.0060844 

29 

3d®(377)4s 

‘^Hll/2 

0.1952878 

4 

3(l®(®D)4s 

®-D3/2 

0.0078607 

30 

3d®(377)4s 

‘^Ps/2 

0.1966664 

5 

3(l®(®D)4s 


0.0089036 

31 

3d®(377)4s 

"‘777/2 

0.1978536 

6 

3d’^ 

^F^/2 

0.0170641 

32 

3d®(3p2)45 

^P9/2 

0.2062854 

7 

Sd’’ 

^^7/2 

0.0221447 

33 

3d®(3p2)45 

"‘-^7/2 

0.2078633 

8 

3d’^ 

^-^5/2 

0.0258613 

34 

3d®(3p2)45 

"‘^’ 6/2 

0.2090388 

9 

3d’^ 

^^3/2 

0.0284084 

35 

3d®(3p2)45 

^P3/2 

0.2098767 

10 

3d®(®D)4s 

'‘P7/2 

0.0724940 

36 

3d®4s2 

®S5/2 

0.2124859 

11 

3d®(®D)4s 

“*^5/2 

0.0764730 

37 

3d®(3G)4s 

^ Gii /2 

0.2317241 

12 

3d®(®D)4s 

“*^3/2 

0.0791021 

38 

3d®(3G)4s 


0.2351555 

13 

3d®(®D)4s 

^ Oi /2 

0.0806177 

39 

3d®(3G)4s 

"‘G'7/2 

0.2367620 

14 

3d’^ 

*P^/2 

0.1227879 

40 

3d®(3G)4s 

"^^5/2 

0.2374345 

15 

3d’^ 

^P3/2 

0.1245992 

41 

3d®(3p2)45 

^P3/2 

0.2349939 

16 

3d’^ 

^Pl/2 

0.1267101 

42 

3d®(3p2)45 

2 p , 
^ 1/2 

0.2454293 

17 

3d’^ 


0.1443871 

43 

3d®(377)4s 

'^Pll/2 

0.2384802 

18 

3d’^ 

'^G-j/2 

0.1491686 

44 

3d®(377)4s 

^^9/2 

0.2401441 

19 

3d’^ 

^P3/2 

0.1673145 

45 

3d®(3p2)45 

"f’7/2 

0.2489119 

20 

3d’^ 

^Pl/2 

0.1721090 

46 

3d®(3p2)45 

= ^5/2 

0.2516957 

21 

3d’^ 

‘^Pll/2 

0.1853545 

47 

3d®(3G)4s 

^^9/2 

0.2769208 

22 

3d’^ 

^ P'S/2 

0.1895961 

48 

3d®(3G)4s 

^G-j/2 

0.2803466 

23 

3d’^ 

"-D5/2 

0.1869643 

49 

3d®(3p)45 

"‘^3/2 

0.2858138 

24 

3d’^ 

^^3/2 

0.1941731 

50 

3d®(3p)45 

"‘O1/2 

0.2858504 

25 

3d®(®P2)4s 

^P3/2 

0.1898222 

51 

3d®(3p)45 

"‘-D5/2 

0.2860280 

26 

3d®(®P2)4s 

^P3/2 

0.1987661 

52 

3d®(3p)45 

^^1/2 

0.2868958 


Note. — * From 


NIST v.4.1 


I 20111'! 
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Table 4. Theoretical radiative widths (s for even-parity levels 


Level 

SST“ 

HFR*’ 

HFR'^ 

CIV3‘* 

TFDAc= 

Q96f 

7-config^ 

TFDAc'^ 

Rec® 

Unc^ 

2 

2.13-3 

2.05-3 

2.14-3 

2.14-3 

2.13-3 

2.13-3 

2.13-3 

2.13-3 

2.12-3 

1.31 

3 

1.57-3 

1.56-3 

1.58-3 

1.58-3 

1.57-3 

1.57-3 

1.57-3 

1.57-3 

1.57-3 

0.38 

4 

7.19-4 

7.28-4 

7.21-4 

7.21-4 

7.18-4 

7.18-4 

7.18-4 

7.18-4 

7.20-4 

0.46 

5 

1.89-4 

1.94-4 

1.93-4 

1.89-4 

1.88-4 

1.88-4 

1.88-4 

1.88-4 

1.89-4 

1.10 

6 

9.99-5 

5.83-5 

6.49-5 

1.42-4 

8.08-6 

1.92-5 

1.62-5 

6.78-5 

6.8 -5 

30 

7 

5.95-3 

6.07-3 

5.94-3 

6.01-3 

5.84-3 

5.85-3 

5.85-3 

5.90-3 

5.92-3 

1.32 

8 

3.99-3 

4.21-3 

3.99-3 

4.03-3 

3.92-3 

3.93-3 

3.92-3 

3.95-3 

3.99-3 

2.29 

9 

1.44-3 

1.55-3 

1.44-3 

1.46-3 

1.41-3 

1.42-3 

1.42-3 

1.43-3 

1.44-3 

2.92 

10 

1.43-2 

1.46-2 

1.43-2 

1.40-2 

1.79-2 

1.77-2 

1.22-2 

1.39-2 

1.48-2 

12.6 

11 

1.38-2 

1.40-2 

1.28-2 

1.30-2 

1.67-2 

1.48-2 

1.19-2 

1.35-2 

1.38-2 

9.90 

12 

1.24-2 

1.27-2 

1.19-2 

1.16-2 

1.55-2 

1.36-2 

1.07-2 

1.23-2 

1.26-2 

10.8 

13 

1.49-2 

1.55-2 

1.08-2 

1.08-2 

1.53-2 

1.31-2 

1.01-2 

1.16-2 

1.26-2 

16.8 

14 

5.65-2 

5.15-2 

4.72-2 

5.24-2 

5.10-2 

4.32-2 

4.13-2 

5.38-2 

5.01-2 

10.2 

15 

5.06-2 

4.73-2 

4.27-2 

4.76-2 

4.86-2 

4.12-2 

3.97-2 

4.53-2 

4.54-2 

7.92 

16 

5.04-2 

4.76-2 

4.27-2 

4.76-2 

4.98-2 

4.15-2 

3.97-2 

4.48-2 

4.54-2 

8.16 

17 

1.94-1 

2.02-1 

1.93-1 

2.52-1 

2.59-1 

1.42-1 

1.36-1 

2.45-1 

2.08-1 

22.4 

18 

1.05-1 

1.11-1 

1.06-1 

1.38-1 

1.39-1 

7.69-2 

7.36-2 

1.32-1 

1.13-1 

22.2 

19 

1.64-1 

1.59-1 

1.60-1 

2.20-1 

2.52-1 

1.19-1 

1.33-1 

2.70-1 

1.94-1 

30.5 

20 

9.58-2 

9.43-2 

1.00-1 

1.27-1 

1.34-1 

7.55-2 

8.55-2 

1.42-1 

1.11-1 

23.1 

21 

1.48-2 

1.58-2 

1.60-2 

1.95-2 

1.20-1 

1.25-2 

1.21-2 

1.45-2 

1.50-2 

15.4 

22 

6.02-2 

6.42-2 

6.44-2 

7.84-2 

3.16-2 

4.69-2 

4.33-2 

9.95-2 

6.07-2 

32.9 

23 

3.78-1 

3.77-1 

4.07-1 

5.30-1 

5.91-1 

2.55-1 

2.35-1 

5.73-1 

4.35-1 

31.6 

24 

5.03-1 

5.15-1 

4.94-1 

6.69-1 

7.53-1 

3.62-1 

3.34-1 

7.25-1 

5.65-1 

28.3 

25 

1.09+0 

1.20+0 

1.11+0 

1.17+0 

1.49+0 

1.15+0 

1.24+0 

9.92-1 

1.16+0 

13.0 

26 

1.30+0 

1.44+0 

1.41+0 

1.45+0 

1.72+0 

1.37+0 

1.43+0 

1.56+0 

1.43+0 

10.6 

27 

1.40+0 

1.56+0 

1.53+0 

1.60+0 

1.87+0 

1.47+0 

1.57+0 

1.28+0 

1.51+0 

12.0 

28 

5.56-1 

6.05-1 

5.40-1 

4.76-1 

6.67-1 

4.77-1 

4.85-1 

5.43-1 

5.45-1 

12.4 

29 

5.21-1 

5.60-1 

5.05-1 

4.47-1 

8.00-1 

4.50-1 

4.58-1 

5.16-1 

5.30-1 

20.4 

30 

4.93-1 

5.25-1 

4.78-1 

4.30-1 

6.37-1 

4.28-1 

4.36-1 

4.86-1 

4.89-1 

13.2 

31 

4.68-1 

4.94-1 

4.52-1 

4.02-1 

1.24+0 

4.06-1 

4.13-1 

4.62-1 

4.45-1 

7.58 

32 

1.03+0 

1.14+0 

1.13+0 

1.21+0 

1.20+0 

1.08+0 

1.05+0 

9.54-1 

1.08+0 

8.79 

33 

8.77-1 

9.76-1 

9.63-1 

1.01+0 

1.03+0 

9.06-1 

8.95-1 

8.22-1 

9.22-1 

8.29 
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Table 4—Continued 


Level 

SST“ 

HFR*> 

HFR'^ 

ClVO-i 

TFDAc= 

Q,96f 

7-configS' 

TFDAc'^ 

Rec® 

Unc.’ 

34 

7.15-1 

7.98-1 

7.81-1 

8.08-1 

8.59-1 

7.17-1 

7.26-1 

6.83-1 

7.52-1 

8.21 

35 

5.76-1 

6.44-1 

6.15-1 

6.21-1 

7.06-1 

5.46-1 

5.76-1 

5.60-1 

6.01-1 

8.54 

36 

3.81+0 

4.54+0 

4.98+0 

4.29+0 

5.57+0 

4.68+0 

4.86+0 

2.30+0 

4.15+0 

27.8 

37 

1.29+0 

1.42+0 

1.29+0 

1.04+0 

1.32+0 

6.73-1 

1.23+0 

7.98-1 

1.10+0 

25.2 

38 

1.32+0 

1.44+0 

1.34+0 

1.16+0 

1.54+0 

7.60-1 

1.24+0 

8.97-1 

1.18+0 

22.9 

39 

1.30+0 

1.41+0 

1.33+0 

1.16+0 

1.61+0 

1.16+0 

1.21+0 

1.29+0 

1.31+0 

10.7 

40 

1.58+0 

1.37+0 

1.30+0 

1.14+0 

1.59+0 

1.15+0 

1.20+0 

1.26+0 

1.32+0 

12.8 

41 

5.14-1 

5.33-1 

5.71-1 

6.83-1 

5.32-1 

6.19-1 

5.32-1 

4.15-1 

5.35-1 

16.1 

42 

6.30-1 

6.67-1 

7.06-1 

8.33-1 

6.44-1 

7.36-1 

6.75-1 

5.17-1 

6.58-1 

15.2 

43 

1.52-1 

1.92-1 

1.96-1 

2.69-1 

1.31-1 

7.12-1 

1.01-1 

7.29-1 

1.74-1 

34.1 

44 

6.61-2 

8.23-2 

8.23-2 

8.65-2 

3.83+0 

5.62-1 

5.99-2 

7.31-1 

7.54-2 

15.5 

45 

5.21-1 

5.41-1 

5.91-1 

6.57-1 

5.90-1 

6.16-1 

5.67-1 

4.32-1 

5.64-1 

12.1 

46 

3.28-1 

3.41-1 

3.75-1 

3.95-1 

3.63-1 

3.90-1 

3.41-1 

2.86-1 

3.45-1 

11.8 

47 

2.98-1 

3.20-1 

3.00-1 

2.86-1 

5.18-1 

3.00-1 

2.81-1 

2.80-1 

2.92-1 

5.19 

48 

2.75-1 

2.92-1 

2.75-1 

2.76-1 

3.49-1 

2.68-1 

2.58-1 

2.82-1 

2.84-1 

9.19 

49 

1.60+0 

1.81+0 

1.73+0 

1.57+0 

2.04+0 

1.57+0 

1.59+0 

1.54+0 

1.67+0 

9.95 

50 

1.62+0 

1.82+0 

1.76+0 

1.60+0 

2.04+0 

1.61+0 

1.61+0 

1.57+0 

1.69+0 

9.34 

51 

1.62+0 

1.83+0 

1.75+0 

1.58+0 

2.05+0 

1.58+0 

1.60+0 

1.59+0 

1.69+0 

9.6 

52 

1.76+0 

1.99+0 

1.91+0 

1.73+0 

2.21+0 

1.71+0 

1.74+0 

1.72+0 

1.83+0 

9.41 

Note. — The level radiative width is 

mended values and their uncertainties, a 

defined as 

± 6 = a X 

Ai = J2j 

io±b 

A{i j). 

The last two columns present our 

recom- 


=‘SST(QDZ96) 
'=HFR(QDZ96) 
'^HFR new 
'*CIV3(DH11) 

®BP extend TFDAc 
^Q96+4d^-corr 


S7-config 







hNewTFDAc 
^Recommended value 


■^Uncertainty (%) 



Table 5. Theoretical branching ratios for forbidden transitions among even parity levels 


i 

j 

SST“ 

HFR*' 

HFR" 

CIV3‘* 

TFDAc" 

Q,96f 

7-config^ 

TFDAc'' 

Rec* 

Unc^ 

2 

1 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

0.00 

3 

2 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

0.00 

4 

3 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

0.00 

5 

4 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

1.00+0 

0.00 

6 

1 

9.16-1 

9.18-1 

9.21-1 

9.15-1 

9.29-1 

9.32-1 

9.32-1 

9.06-1 

9.21-1 

1.01 

6 

2 

8.37-2 

8.18-2 

7.80-2 

8.17-2 

7.00-2 

6.72-2 

6.79-2 

9.40-2 

7.80-2 

11.9 

7 

6 

9.82-1 

9.90-1 

9.88-1 

9.75-1 

9.98-1 

9.97-1 

9.97-1 

9.88-1 

9.89-1 

0.82 

8 

7 

9.82-1 

9.90-1 

9.87-1 

9.75-1 

9.97-1 

9.95-1 

9.97-1 

9.90-1 

9.89-1 

0.78 

9 

8 

9.79-1 

9.87-1 

9.86-1 

9.73-1 

1.00+0 

9.93-1 

9.93-1 

9.86-1 

9.87-1 

0.86 

10 

1 

3.31-1 

3.55-1 

3.95-1 

3.76-1 

2.55-1 

4.62-1 

3.59-1 

3.49-1 

3.60-1 

16.3 

10 

2 

9.16-2 

9.86-2 

1.10-1 

1.06-1 

7.71-2 

1.36-1 

1.12-1 

1.03-1 

1.04-1 

16.4 

10 

3 

5.89-2 

6.12-2 

6.88-2 

6.94-2 

4.68-2 

8.64-2 

6.78-2 

6.15-2 

6.51-2 

17.4 

10 

6 

4.18-1 

3.92-1 

3.44-1 

3.62-1 

4.84-1 

2.53-1 

3.72-1 

3.90-1 

3.77-1 

17.5 

10 

7 

9.23-2 

8.56-2 

7.55-2 

8.00-2 

1.20-1 

5.47-2 

8.28-2 

8.63-2 

8.47-2 

21.5 


Note. — The branching ratio is given by b(i —>■ j) = A(i — >■ j) X T(i) = A{i j)/ A{i — )■ k). The last two 

columns tabulate our recommended values and their uncertainties, a ± b = a X 10^*”. 

'*SST(QDZ96) 

‘’HFR(QDZ96) 

‘^HFR new 

‘^CIV3(DH11) 

®BP extend TFDAc 

^Q96+4d^-corr 

®7-config 

•'NewTFDAc 

'Recommended value 


j Uncertainty (%) 
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Table 6. 


Comparison of theoretical lifetimes (s with experiment 


Level 

SST“ 

HFR*’ 

HFR'^ 

CIV3‘* 

TFDAc® 

Q96^ 

7-config^ 

TFDAc'^ 

Rec* 

Expt-^ 

36 

0.262 

0.220 

0.201 

0.233 

0.180 

0.214 

0.206 

0.435 

0.24 ±0.07 

0.23 ± 0.03 

37 

0.775 

0.704 

0.775 

0.962 

0.758 

1.49 

0.813 

1.25 

0.91 ±0.23 

0.75 ±0.01 

38 

0.758 

0.694 

0.746 

0.862 

0.649 

1.32 

0.806 

1.11 

0.85 ±0.20 

0.65 ± 0.02 

43 

6.58 

5.21 

5.10 

3.72 

7.63 

1.40 

9.90 

1.37 

5.8 ±2.0 

3.8 ±0.3 


=='SST(QDZ96) 
‘>HFR(QDZ96) 
'^HFR new 


'iCIV3(DHll) 

®BP extend TFDAc 
fQ96+4d^-corr 
®7-config 
•'NewTFDAc 


‘Recommended theoretical lifetimes 


JExperimental lifetimes from the Ferrum Project ( Hartman et al. 


2003 


Gurell et al. 


2009f 
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Table 7. Comparison of theoretical and observed branching ratios 


Calculation 


SST(QDZ96) 

5.83 

HFR(QDZ96) 

5.09 

HFR new 

5.16 

CIV3(DH11) 

4.76 

BP extend TFDAc 

358 

Q96+4d^-corr 

12.64 

7-config 

4.10 

NewTFDAc 

4.92 

Recom.({5/?th = 0)^ 

4.34 

Recommended^ 

1.38 


Note. — The comparison (in 
terms of the reduced x^) is for 
106 observed branching ratios 
with common upper levels. 

^Excludes estimated theoret¬ 
ical uncertainties. 

'^Includes estimated theoret¬ 


ical uncertainties. 
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Table 8. Effective collision strengths from the ground level of [Fe II] at 10^ K 


j 

Q+RM“ 

Q+RM*' 

Q+RM'^ 

Q+RM'^ 

7-config^ 

DARC-'' 

Mean9 

o-'* 

ZP96* 

BP96'= 

RH07'' 

2 

1.81+0 

2.05+0 

2.80+0 

2.24+0 

2.85+0 

5.16+0 

2.31+0 

16.0 

5.52+0 

4.65+0 

4.84+0 

3 

3.36-1 

3.89-1 

6.52-1 

4.93-1 

4.20-1 

1.20+0 

4.21-1 

23.8 

1.49+0 

1.29+0 

1.12+0 

4 

1.54-1 

1.58-1 

3.58-1 

2.33-1 

1.78-1 

5.11-1 

1.95-1 

35.3 

6.84-1 

8.13-1 

5.29-1 

5 

7.16-2 

7.33-2 

1.81-1 

1.16-1 

8.57-2 

2.22-1 

9.51-2 

37.8 

2.84-1 

4.33-1 

2.47-1 

6 

2.05+0 

1.99+0 

2.56+0 

2.20+0 

1.84+0 

3.07+0 

2.05+0 

11.5 

3.60+0 

1.31+0 

2.83+0 

7 

7.46-1 

7.53-1 

9.27-1 

8.74-1 

6.78-1 

1.33+0 

7.64-1 

11.6 

1.51+0 

6.14-1 

1.21+0 

8 

1.57-1 

1.59-1 

1.90-1 

2.26-1 

1.42-1 

4.08-1 

1.64-1 

18.1 

4.97-1 

1.35-1 

2.98-1 

9 

3.18-2 

2.72-2 

3.29-2 

6.63-2 

2.25-2 

1.01-2 

3.09-2 

48.1 

1.37-1 

1.65-2 

5.84-2 

10 

1.61+0 

1.59+0 

2.02+0 

1.78+0 

4.26+0 

9.04+0 

2.37+0 

45.5 

1.10+1 

1.43+1 

1.04+1 

11 

1.18-1 

1.02-1 

1.47-1 

1.51-1 

1.17-1 

6.22-1 

1.11-1 

24.2 

5.59-1 

5.72-1 

5.26-1 

12 

5.63-2 

4.45-2 

6.92-2 

8.79-2 

4.27-2 

2.14-1 

5.13-2 

34.1 

1.91-1 

3.86-1 

1.79-1 

13 

2.21-2 

2.19-2 

3.01-2 

3.89-2 

1.94-2 

7.61-2 

2.33-2 

31.0 

6.01-2 

2.33-2 

6.49-2 

14 

7.53-1 

6.40-1 

8.90-1 

7.30-1 

6.75-1 

9.42-1 

7.02-1 

11.0 

9.48-1 

5.42-1 

9.47-1 

15 

3.69-1 

3.78-1 

4.78-1 

3.91-1 

4.06-1 

4.80-1 

3.98-1 

8.3 

5.02-1 

3.19-1 

5.22-1 

16 

1.19-2 

7.60-3 

9.28-3 

9.40-3 

5.46-3 

5.89-2 

7.44-3 

21.8 

3.28-2 

8.20-3 

6.11-2 

17 

2.35-1 

7.73-2 

1.72-1 

2.01-1 

3.42-2 

9.79-2 

1.08-1 

59.1 

- 

- 

9.73-2 

18 

1.24-1 

4.35-2 

6.24-2 

1.12-1 

1.28-2 

5.65-2 

5.19-2 

65.4 

- 

- 

6.25-2 

19 

2.84-2 

1.78-2 

2.76-2 

3.04-2 

8.05-3 

2.98-2 

1.77-2 

55.5 

- 

- 

4.54-2 

20 

4.52-3 

7.12-5 

3.22-3 

6.95-3 

3.70-4 

8.04-3 

1.84-2 

244. 

- 

- 

1.32-2 

21 

8.65-2 

3.39-5 

1.18-1 

1.02-1 

4.08-3 

4.46-1 

4.23-2 

108. 

- 

- 

1.55-1 

22 

6.33-2 

1.53-3 

5.04-2 

7.77-2 

1.79-3 

1.74-1 

2.59-2 

107. 

- 

- 

3.90-2 

23 

3.86-2 

3.02-1 

3.12-2 

5.29-2 

3.26-3 

4.60-2 

8.80-2 

114. 

- 

- 

4.66-2 

24 

9.19-3 

9.82-2 

6.72-3 

1.50-2 

8.14-4 

9.25-3 

2.27-2 

1430 

- 

- 

1.63-2 

25 

2.30-1 

1.23-1 

3.35-1 

2.85-1 

1.44-1 

3.89-1 

1.85-1 

56.0 

3.08-1 

2.15-1 

2.67-1 

26 

5.85-2 

4.14-2 

8.79-2 

6.86-2 

3.49-2 

1.02-1 

6.48-2 

53.7 

9.92-2 

5.17-2 

7.01-2 

27 

1.78-3 

2.26-4 

2.62-3 

2.14-3 

9.46-4 

1.69-2 

3.67-3 

172. 

8.00-3 

4.30-3 

6.12-3 

28 

5.42-1 

5.66-1 

8.14-1 

6.88-1 

3.26-1 

7.97-1 

5.44-1 

30.2 

6.31-1 

1.96-1 

7.19-1 

29 

3.13-1 

3.45-1 

4.57-1 

4.50-1 

1.77-1 

2.70-2 

3.14-1 

33.6 

3.11-1 

9.55-2 

3.74-1 

30 

1.18-1 

1.39-1 

1.71-1 

1.88-1 

6.60-2 

6.76-1 

1.19-1 

37.5 

9.51-2 

2.76-2 

1.46-1 
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Note. — a ±b = a X 10=*=^. 

'"Q+RM-ns 

*’Q+RM-shift 

‘'Q+RM-ns-RA=8 

‘iQ+RM-ns-RA=14.5 

®7-config 

fDARC 


®Mean of effective collision strengths obtained in the present calculations 

^Standard deviation (%) for effective collision strengths obtained in the present calculations 


Zhane and Pradhan 

1199S1 


Bautista Sz Pradhan 

11996) 


^ Ramsbottom et al 


I 20071^ 
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Table 9. Comparison of normalized line intensities in HH 202 with model predictions 


Model 


Te 

(10^ K) 

rie 

(10^ cm~^) 

ZP96“ 

1.40 

0.9 

4.7 

BP98*’ 

1.39 

1.0 

4.5 

RH07= 

1.60 

1.2 

3.0 

Pres"^ 

3.13 

1.0 

2.0 

Pres*" 

1.28 

1.2 

10. 

Pres-^ 

1.05 

1.0 

7.0 

Pres® 

1.01 

1.15 

6.6 


Note. — Each model uses a specific set 
of effective collision strengths. The listed re¬ 
duced corresponds to the best match with 
observation at the quoted temperature and 
density. 


^Collision 


strengths 

Zhane and Pradhan 

(1995 

). 

^Collision 


strengths 

Bautista Sz Pradhan 

1998 

). 

‘^Collision 


streng 

ths 

Ramsbottom et al. 



2007) 

Bautista Sz Pradhan 

1998 

). 


^Present DARC calculation. 


from 


from 


from 

and 


^Present bprm calculation with 114 levels. 
^Present bprm calculation with 63 levels. 


^Present bprm calculation with 7-config. 

























